Method for operation of hydropower reservoir with a 2-parameter elevation rule curve

ABSTRACT

Disclosed herein is a novel method for operating a hydropower reservoir which is an improvement over the existing single-parameter (the current month) USACE Rule Curve approach, the improvement comprising the consideration of a second parameter, namely the water level of the reservoir at the beginning of the month, in the decision-making process for operation of the reservoir.

CLAIM OF PRIORITY

This application is a continuation of U.S. patent application Ser. No. 15/408,855 filed Jan. 18, 2017, now U.S. Pat. No. 10,125,740, which claims the benefit of U.S. Provisional Patent Application Ser. No. 62/279,943 filed on Jan. 18, 2016, and U.S. Provisional Patent Application Ser. No. 62/331,648 filed on May 4, 2016, the contents of which are incorporated herein by reference.

FIELD OF INVENTION

The present invention relates to optimizing the operation of and power collected from a hydro powered reservoir.

BACKGROUND OF THE INVENTION

Many studies on planning and operation of reservoirs for power production rely on the operational strategy from the United States Army Corps of Engineers (“USACE”) described in reference EM-1110-2-1701 entitled “Engineering and Design, Hydropower” (the “USACE Reference”), which aims to preserve the firm energy of the reservoir and then to produce secondary energy. The operational strategy presented is based on two aspects: (a) the actual control strategy, and (b) a set of control values of the reservoir level, which acts as a threshold in the decision-making process, and is called the Rule Curve. This Rule Curve has one parameter, which is the current month of the year (it is a 1-parameter set), and gives a control value of the reservoir level at the end of each calendar month.

The operational strategy to use with this set of control parameters is simple, stated as follows: (i) at the beginning of the current month simulate the operation of the reservoir/powerhouse generating only the continuous power and firm energy; (ii) if the simulated end-of-month level in the reservoir falls below the control value, generate only the firm energy and continuous power; (iii) if the simulated end-of-month level is greater than the control value, generate as much energy as possible until the level reaches down to the control value, always limited by the power-plant capacity. The set of control values is fixed over the years and gives one control value of the end-of-month reservoir level for each calendar month.

The Rule Curve and the continuous power/firm energy yield are obtained by a recursive backward simulation algorithm until the maximum firm yield and continuous power are obtained. The calculation of the continuous power/firm energy yield is made by trial-and-error, generating only the trial continuous power under the assumed daily operation pattern, and calculating the required level at the beginning of each time period considering the reservoir at the minimum power pool level at the end of the time series; the reservoir fails if it requires a level greater than the maximum power pool level. This backward simulation is repeated until the maximum continuous power/firm energy yield is obtained.

The Rule Curve is the envelope of the maximum end-of-month levels obtained in the backward simulation which gives the firm energy yield and continuous power. As a by-product of this procedure, at least one value of the Rule Curve is the top of the power pool. The Rule Curve can be updated periodically as more hydrological information is collected during the actual operation of the reservoir, or bathymetric studies are undertaken to analyze sediment deposition and update the reservoir elevation-area-storage curve, or other significant event occurs.

FIG. 1 shows a typical Rule Curve (RC) for a reservoir, following the USACE guidelines. It shows also the minimum and maximum power pool levels, and the general recommendation for operation of the power-plant. Note that one level of the Rule Curve equals the maximum normal power pool level. As can be seen, under the USACE operational strategy, the control of the reservoir operation depends only on the current time of the year. There is no influence of the current state of the reservoir system in the selection of the value of the end-of-month control level, nor is there an explicit maximization of the mean energy generated by the reservoir.

SUMMARY OF THE INVENTION

Disclosed herein is a new approach to generate an improvement for the Rule Curve, to be used with the control strategy already stated in the USACE Reference. The improvement is the inclusion of a parameter which represents the actual state of the reservoir system in the decision-making process. The parameters used, then, are current time of the year and current reservoir level (at the beginning of the month). This approach is a two-tiered scheme:

1) The estimation of the value of a composite parameter of the reservoir operation, ΔV H_(ef) (by month or any suitable time-period), using the available hydrological data (inflow, rainfall, evaporation), for the generation of continuous power and firm energy; this composite parameter represents the monthly maximum effective potential energy withdrawn from (negative values) or stored into the reservoir (positive values). The calculations are done with additional data processing at each time step to calculate the monthly variation of the effective potential energy withdrawn from or stored into the reservoir. At the end of the simulation, the algorithm takes the upper envelope of the maximum monthly values of the composite potential energy parameter, instead of the envelope of the maximum values of the end-of-month storage.

2) The use of the composite monthly maximum parameter ΔV H_(ef) to estimate the end-of-month control reservoir levels. This task is done by solving the mass balance and net head equations with the appropriate efficiency, to give an end-of-month level to any given beginning-of-month reservoir level within the range of admissible levels in the reservoir power pool. Since the mass-balance equations involve the effective rainfall and evaporation over the reservoir surface, the algorithm takes a conservative approach and evaluates the mass balance with the minimum monthly values of the difference between the effective rainfall and effective evaporation. This procedure provides a one-to-one relationship between the initial reservoir level and the recommended final level for any month. In real-time operations, this calculation can be made with the best estimate of the rainfall and evaporation over the reservoir, to produce the end-of-month Rule Curve levels, to further enhance the estimate of the end-of-month reservoir levels.

An example of the 2-Parameter Elevation Rule Curve is presented in FIG. 2, in which the end-of-month level is given as a function of the beginning-of-month level and the time of the year and the minimum monthly values of the difference between effective rainfall and evaporation over the reservoir surface. The rule curve ranges from the minimum to the maximum reservoir power pool levels. For clarity purposes, it is presented in four seasons, each with duration of three months.

The calculations are simple and straightforward in the derivation of the composite parameter and the estimation of the Rule Curve. There is no maximization algorithm, only the consideration of a parameter representative of the current state of the system. The use of this 2-Parameter Rule Curve with the control strategy presented by USACE can increase the annual energy output of the reservoirs while maintaining the firm energy yield. The increase in total energy can range from 2-5% up to 20-25% and 50-60% in some applications, depending on the characteristics of the reservoir and power-plant system, such as storage, rated capacity, net head variations and inflow.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 depicts a typical USACE hydropower Rule Curve.

FIG. 2 depicts a typical example of the 2-Parameter Elevation Rule Curve.

FIG. 3 depicts the operating pattern used to simulate the generation of continuous power and firm energy.

FIG. 4 depicts a tabulation of the hydrological parameters of the Tayucai Reservoir.

FIG. 5 depicts a tabulation of the Reservoir Level-Area-Storage Curve for the Tayucai Reservoir.

FIG. 6 depicts a tabulation of the Tailwater Discharge Levels for the Tayucai Power-Plant.

FIG. 7 depicts a tabulation of the Continuous Power Output and Finn Energy Yield for a Backward Simulation at the Tayucai Reservoir, Design Head 80 in.

FIG. 8 depicts a tabulation of the Maximum Parameters PTX at the Tayucai Reservoir; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 9 depicts the estimation of the 2-Parameter Rule Curve for the simulation of Unrestricted A10 Case 1 for Tayucai Reservoir; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 10 depicts the estimation of the 2-Parameter Rule Curve for this simulation of Restricted A20 Case 1 for Tayucai Reservoir; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 11 depicts the estimation of the 2-Parameter Rule Curve for this simulation of Unrestricted A10 Case 2 for Tayucai Reservoir; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 12 depicts the estimation of the 2-Parameter Rule Curve for this simulation of Restricted A20 Case 2 for Tayucai Reservoir; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 13 depicts a tabulation of the Tayucai Reservoir USACE Rule Curve, End-of-Month Storage and Reservoir Levels; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 14 depicts the USACE hydropower Rule Curves for Tayucai reservoir, Cases 1 and 2; Installed Capacity 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 15 depicts the 2-Parameter Unrestricted A10 and USACE Rule Curves for Tayucai Reservoir: 3,800 MW, Case 1, Month of March; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 16 depicts a combined 2-Parameter and USACE Rule Curves Unrestricted A10 for Tayucai Reservoir; 3,800 MW, Case 1, Month of March; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 17 depicts several linear combinations of the Unrestricted A10 Rule and the USACE Rule Curves, on level and storage with weights 0.25, 0.5, and 0.75, as well as the original A10 and USACE Rule Curves; 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 18 depicts a tabulation of Mean Annual Energy, Combined Extreme Values, 2-Parameter and USACE Rule Curves for Tayucai Reservoir; 3,800 MW; Design Head 80 mca; Design Flow 5,667 m3/s.

FIG. 19 depicts a tabulation of Mean Annual Energy, Linear Combination of Levels, 2-Parameter and USACE Rule Curves for Tayucai Reservoir, 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 20 depicts a tabulation of Mean Annual Energy, Linear Combination of Storage, 2-Parameter and USACE Rule Curves for Tayucai Reservoir, 3,800 MW; Design Head 80 in; Design Flow 5,667 m3/s.

FIG. 21 depicts a tabulation of the Tayucai Powerplant Characteristics, for diverse Installed Capacities, Design Head and Design Flow.

FIG. 22 depicts a tabulation of the Tayucai Reservoir Energy Generated by the use of the USACE Rule Curve.

FIG. 23 depicts a tabulation of the Tayucai Reservoir Energy Generated by the Use of the 2-Parameter Unrestricted A10 Rule Curve and Combinations.

FIG. 24 depicts a tabulation of the Tayucai Reservoir Energy Generated by the Use of the 2-Parameter Restricted A20 Rule Curve and Combinations.

FIG. 25 depicts a tabulation of the Tayucai Reservoir Comparison of Energy Generated, by the Use of the USACE and 2-Parameter Rule Curves.

FIG. 26 depicts Mean Annual Energy Generated by use of the Diverse Rule Curves at Tayucai Reservoir.

FIG. 27 depicts Secondary Annual Energy Generated by use of the Diverse Rule Curves at Tayucai Reservoir.

DETAILED DESCRIPTION OF INVENTION

Analytical Description of the Reservoir Power Plant System

Operation Pattern for Continuous Power/Firm Energy

FIG. 3 shows the operating pattern used to simulate the generation of continuous power and firm energy. The pattern follows a base/peak load, with duration Tb and Tp, respectively. The Continuous Power is indexed to the base/peak period, CPb and CPp. The restrictions on the Continuous Power and pattern are the following:

a) CPp and CPb must be lower or equal than the power output at the minimum power pool level.

b) Total time, Tb+Tp, equals 24 hr.

c) The relationship between the base and peak continuous power, fB=CPb/CPp, lies between zero and one. A value of zero means the power-plant is dedicated to generate peak load, whereas a value of one means the plant generates in base load only.

Mass Balance of Reservoir

The mass balance of the reservoir is represented by the following recursive equations, for any given period: Vf=Vi+ΔV  (“EQ 1”)

$\begin{matrix} {{\Delta\; V} = {{{Qinf}\;\Delta\; T} - {{Qtp}\;\Delta\;{Tp}} - {{Qtb}\;\Delta\;{Tb}} + {\left( {{{Kr}\; R} - {{Ke}\;{Ev}}} \right)\;{Area}} - {{Qeco}\;\Delta\; T} - {{Qloss}\;\Delta\; T}}} & \left( {``{{EQ}\mspace{14mu} 2}"} \right) \end{matrix}$

-   -   Subject to the inequality constrain:         Vmin≤Vf≤Vmax  (“EQ 3”)

In which:

-   -   Vf: Volume stored in reservoir at the end of the month.     -   Vi: Volume stored in reservoir at the beginning of the month,         known value.     -   ΔV: Variation on volume stored during the period, without         restrictions imposed by the reservoir levels (maximum and         minimum).     -   Qinf: Average Inflow into the reservoir during the current         period.     -   R: Total monthly rainfall over the reservoir area during the         current period.     -   Ev: Total monthly evaporation over the reservoir area during the         current period.     -   Qtp: Average monthly flow into turbines during peak time, Tp.     -   Qtb: Average monthly flow into turbines during base time, Tb.     -   Qeco: Average ecological flow and other demands during the         period which do not pass through the power house.     -   Qloss: Average water losses during the period.     -   Vmin: Volume stored at the minimum reservoir power pool level.     -   Vmax: Volume stored at the maximum reservoir power pool level.     -   ΔT: Time conversion factor from the units of flow into total         monthly volume, for full day.     -   ΔTp: Time conversion factor from the units of flow into total         monthly volume, for peak hours.     -   ΔTb: Time conversion factor from the units of flow into total         monthly volume, for base hours.     -   Area: Average surface area of the reservoir during the period.     -   Kr: Effective rainfall and conversion constant.     -   Ke: Effective evaporation and conversion constant.

The relationship between level-area storage of the reservoir is represented by the level-area-storage curve. If the water surface level falls below the minimum power pool level, the reservoir cannot yield the required power and energy output, and the available water must be allocated to the diverse demands on a predefined priority basis. The final level is then set to the minimum power pool level. If the water level goes above the maximum power pool level, it generates a spill, which is then calculated. The final level is set to the maximum power pool level. This procedure is repeated on a month-by-month basis, for the duration of the hydrological time series available.

Power Output of Power Plant

The equations representing the power output of the power-plant are the following: P=γQt H ef _(t) ef _(g)  (“EQ 4”) E=PpΔTp+PbΔTb  (“EQ 5”)

-   -   Subject to the inequality constrains on the continuous power to         generate firm energy:         Pp≥CPp  (EQ 6) and         Pb≥CPb  (“EQ 7”)

In which:

-   -   P: Power output of the power-plant, it is indexed for the         generation on peak time, Pp; and base time, Pb.     -   CPp: Continuous power during peak time.     -   CPb: Continuous power during base time.     -   E: Total energy generated during period.     -   γ: specific weight of water.     -   Qt: water flow into the turbines, indexed for peak and base         time, Qtp and Qtb.     -   H: average monthly net head on the turbines, indexed for peak         and base time, Hp and Hb.     -   ef_(t): average efficiency of the turbine units, indexed for         generation at peak and base time.     -   ef_(g): average combined efficiency of the generator and         transformer units, it is presumed constant throughout the         simulation.

Net Head

The net head is a function of the reservoir elevation, the energy losses through the conveyance systems (free-surface and pressure flow) from the reservoir into the turbines, and the tailwater elevation rating, as follows: H=Elev−HL−TW  (“EQ 8”)

In which:

-   -   H: Average net head on the turbines, indexed for peak and base         power.     -   Elev: Average surface level of reservoir during period.     -   HL: Average energy losses during the period, due to conveyance         from reservoir into the power-plant, indexed for peak and base         power. It is a function of total flow into the power-plant and         unit flow into the turbines.     -   TW: Average tailwater level during the period, indexed for peak         and base power, it is a function of total flow into the power         plant, and must take into account any spillway discharges that         affect the power house discharge channel levels.

Efficiency

The efficiency of the turbines is a function of the net head and unit flow on the turbines. The relationship between efficiency-head-flow is part of the data set and can be read from the efficiency hill-diagrams of typical turbines or of the installed turbines. The combined efficiency of the generator and transformer units is generally presumed constant.

Estimation of the Composite Monthly Reservoir Parameter, PT

The backward simulation algorithm used to define the continuous power and firm energy can be represented by the equations: Vi=Vf−ΔV  (“EQ 9”) Pp=CPp  (“EQ 10”) Pb=CPb  (“EQ 11”)

-   -   Subject to the inequality constrain:         Vmin≤Vi≤Vmax  (“EQ 12”)

In which the variables are defined in equations (EQ 1) and (EQ 2), with the initial condition of reservoir at minimum power pool level at the end of the time series. The equations are solved recursively backwards, beginning at the end of the time series data, up to the initial period.

Once the continuous power and firm energy are identified by the backward simulation procedure, the following parameter is calculated for the time-series: PT(i,j)=ΔV(i,j)Hef(i,j)  (“EQ 13”)

In which:

-   -   i, j: indexes for the year and month of the time series.     -   PT(i,j): monthly value of the variation on stored effective         potential energy in the reservoir.     -   ΔV(i,j): monthly value of the variation on storage during the         period.     -   Hef(i,j): monthly average product of net head and efficiency of         the power-plant (turbines, generator and transformers combined).

Estimation of the Variation in Stored Volume of Reservoir

The calculation of ΔV(i,j) involves two variants, defined as: unrestricted variation (A10) and restricted variation (A20) of stored volume, as follows:

a) A10—Unrestricted variation of stored volume in reservoir. It is calculated prior to the imposition of the inequality constrains on the stored volume from (EQ 2), restated below for ease of reference:

Δ V = Qinf Δ T − Qtp Δ Tp − Qtb Δ Tb + (Kr R − Ke Ev) Area − Qeco Δ T − Qloss Δ T.

b) A20—Restricted variation of stored volume in reservoir. It is calculated after the imposition of the inequality constrains on the stored volume (EQ 12): Vmin≤Vi≤Vmax  (“EQ 12”) ΔV=Vf−Vi  (“EQ 14”)

If the inequality constrains are not binding, the two calculated values, unrestricted (A10) and restricted (A20), are equal.

Estimation of the Average Net Effective Head

The net effective head is the product of net head and efficiency of the power-plant. The net head is defined in equation (EQ 8) restated below, and the efficiency is a function of the net head and unit flow in the turbines. HEF=Hef _(t) ef _(g)  (“EQ 15”) H=Elev−HL−TW  (“EQ 8”)

In which:

-   -   HEF: Net effective head, indexed for peak and base power.     -   ef_(t): Efficiency of the turbines, indexed for peak and base         power.     -   ef_(g): Combined efficiency of the generators and transformers,         it is presumed constant throughout simulation.

The average net effective head is a composite measure over the time period, weighted on the duration of peak and base load. The equation is the following: HEFavg=(HEFpTp+HEFbTb)/24  (“EQ 16”)

In which:

-   -   HEFavg: weighted average net effective head over the period.     -   HEFp: net effective head on peak power.     -   HEFb: net effective head on base power.     -   Tp: duration of peak power.     -   Tb: duration of base power.

Data Processing of the Composite Monthly Reservoir Parameter PT to Obtain the Maximum Monthly Values

The control value for the estimation of the Rule Curve is the maximum envelope of the monthly values of the PT-Parameter, that is: PTX(j)=max(PT(ij)),i=1,n  (“EQ 17”)

In which:

-   -   i,j: indexes for the year and month of the time series.     -   PTX(j): maximum value of PT parameter for month j, over the         entire series.     -   n: duration of the time series, years.

This process is done for the two variants of the PT parameter, the unrestricted and restricted case. Each case gives an envelope to be used in the estimation of the 2-Parameter Rule Curve

Estimation of the End-of-Month Level for the 2-Parameter Elevation Rule Curve

The estimation of the end-of-month level proceeds after the extreme maximum monthly values of the PT parameter is calculated. The calculations follow the inverse procedure to that which gave origin to the PT parameter, and are described below:

For any given month or time period:

a) Define initial reservoir level: Any value ranging between the minimum and maximum power pool levels. To obtain the full range of end-of-month storage, the calculations must be repeated with initial level varying between the extremes, with the adequate variation to obtain a good definition of the curve.

b) Define the extreme values of the effective rainfall and evaporation: The parameter RE(i,j)=+Kr R(i,j)−Ke Ev(i,j) is the effective contribution of the rainfall and evaporation to the mass balance equation. Use the time series available to estimate REmin(j), the minimum monthly value of this parameter to use in the mass-balance equations.

c) Estimate the variation in storage to calculate the required end-of-month storage for all the initial reservoir levels required, considering the generation of continuous power and firm energy only. The equations used are, for any given month and initial reservoir level: ΔV=PTX/Hef  (“EQ 18”) Vf=Vi+ΔV−REmin Area+QecoΔT+Qloss ΔT  (“EQ 19”) Vmin≤Vf≤Vmax  (EQ 3)

In which

-   -   ΔV: variation in storage from the beginning to end of month.     -   PTX: monthly extreme value of the PT parameter for the month.     -   Hef: average effective net head for the period, considering the         generation of continuous power and firm energy only.     -   Vi: reservoir storage at the beginning of the month, associated         with the initial level.     -   Vf: reservoir storage at the end of the month.     -   REmin: minimum value of the effective contribution of rainfall         and evaporation, for the month.     -   Area: average surface area of the reservoir during the period.

d) Use the level-area-storage curve to determine the end-of-month level for the Rule Curve. Once the storage at the end of the month is calculated, use the level-area-storage curve to estimate the reservoir level that defines the Elevation Rule Curve.

Sample Calculations

Provided here is the complete calculation of the 2-Parameter Rule Curve, in line with the calculations proposed by the USACE Reference to develop their Rule Curve, are discussed. First, there are the data requirements, then the processing of the data with the backward simulation to generate the control parameters, and, finally, the estimation of the values of the 2-Parameter Rule Curve for each period and range of values within the admissible levels in the reservoir. The experimental site is the Tayucai reservoir, a site under study in the Upper Caroni Basin, in Venezuela.

Data Requirements

For optimal results, the data should be be the most complete and up-to-date available, covering the hydrological aspects, the characteristics of the reservoir and power-plant, and other water demands to and losses from the reservoir.

Hydrology:

The hydrological data includes the total monthly streamflow into the reservoir, as well as the data for direct monthly rainfall to and evaporation from the reservoir. The length of the record period should be as long as possible. The data is taken from the appropriate sources or estimated from the hydrological modeling of the basin. FIG. 4 shows the monthly and annual means of total inflow, and the monthly rainfall and evaporation with annual total. The value of the parameter Kr, effective rainfall is 0.95; the value of the Ke parameter, effective evaporation, is 0.80.

Reservoir:

The data required for the reservoir are the following:

a) Level-area-storage curve: It is taken from the appropriate topographical maps, or from bathymetric studies of existing reservoirs. The data can be presented in tabular and graphical forms. The interpolation between data points can be made with a line segment, or with the aid of a curve-fitting procedure. The level-area-storage curve is presented in FIG. 5.

b) Representative Levels: It is required to determine the maximum and minimum levels of the power pool, EL_(max) and EL_(min), to define the usable storage for regulation purposes. These levels define the maximum and minimum storage, Vmax and Vmin, and the usable storage as the difference between the maximum and minimum storage values. For this sample, the minimum power pool level is 340 masl, the maximum power pool level is 380 masl. Minimum storage is 1,456 mill. m³, and maximum storage 9,419 mill. m³, the usable storage is 7,963 mill. m³, which is roughly 10% of the mean annual inflow to the reservoir.

Power Plant:

a) Operating pattern: The operating pattern is fully defined with the following parameters:

-   -   Duration of peak time, Tp: 4 hour     -   Duration of base time, Tb: 20 hour     -   Ratio between base and peak continuous power, fB: Case 1: 1.00         (base load), Case 2: 0.40

b) Tailwater level: For reaction turbines, the tailwater level is the rating curve (flow vs. level) of the discharge channel. For impulse turbines, the tailwater level is the fixed elevation of the discharge jets into the turbines. In this example, the rating curve of the discharge channel is shown in FIG. 6.

c) Turbines: The turbines are defined by their number, type, design head and flow, the usable range of net head and flow, the ratio of maximum power over rated power, and by the efficiency hill-diagram. The turbines can be analyzed jointly as a single unit or separated, in which case the number of operating turbines varies with the net head, flow and power produced by the power-plant. In any case, the operation of the turbines should be the most efficient.

The power-plant is located at the foot of the dam, with one penstock for each turbine. In this example, design head is fixed at 80 in. The machines are Francis turbines, with the efficiency hill diagrams for a machine of specific speed 344 (metric). Total power output is about 3,800 MW, with several turbines (in the range of 10 to 15 machines). Hence, it is expected that the most efficient operation of the power-plant will provide that the unit flow of the machines will be greater than or nearly equal to the maximum efficiency flow for any given head and total power output. The combined efficiency of the generator and transformer units is presumed constant, equal to 0.96.

d) Head Losses: The power-plant is located at the foot of the dam. The head losses in the penstock are considered constant, equal to 1.5 in.

Other Demands

Other demands include ecological flow (minimum required flow, fisheries, navigation, other uses), and other water demands such as irrigation or urban water supply. In this case, the only demand is maintaining a minimum flow downstream from the power-plant equal to 10% of mean annual inflow, some 250 m3/s, to maintain ecological flow. If the total turbined flow at peak or base times is greater than the minimum required flow of 250 m3/s, no action is taken; however if the turbined flow is smaller than the minimum required, the difference between the minimum required flow and turbined flow is released into the river channel without going through the power plant.

Estimation of the Composite Parameter PT and Maximum Values PTX

Presented here is the backward simulation algorithm, and the calculation of the required parameters, along with the estimation of the USACE Rule Curve, for comparison and combination purposes, as will be seen below. The calculations are made for the two (2) cases indicated: Case 1 (base load) with fB equal to 1.0; and Case 2 (peak load) with fB equal to 0.4.

Backward Simulation:

The result of the backward simulation is presented in FIG. 7, for Case 1 and Case 2, with the continuous power (peak and base) and the firm energy yield estimated for each case.

Composite Parameter PT and Maximum Parameter Value PTX:

FIG. 8 shows the values of the PTX parameters, arranged by month, from the backward simulation. The table shows the maximum value of the parameter, PTX, for Case 1 and Case 2, in the unrestricted (A10) and restricted (A20) variants.

Estimation of the 2-Parameter Rule Curves

FIGS. 9-12 show the estimation of the 2-Parameter Rule Curve for this simulation for both the Unrestricted A10 and the Restricted A20 for Case 1 and Case 2.

Estimation of the USACE Curve

FIG. 13 shows the maximum end-of-month storage and associated reservoir level of the backward simulation, arranged by month, taken from the backward simulation shown for Case 1 and Case 2. In FIG. 14 the maximum values of the end-of-month reservoir levels are reproduced, to form the USACE Rule Curve. As can be seen, the USACE Rule Curves for Cases 1 and 2 are very similar.

Combined Use of the 2-Parameter and USACE Rule Curves

The principal purpose of both the 2-Parameter and USACE Rule Curves is to preserve the Finn Energy yield of the reservoir. Both elevation rule curves represent a recommended end-of-month level based on the current time of year and, for the 2-Parameter Rule Curve, on the reservoir level at the beginning of the month. FIG. 15 shows a typical 2-Parameter Rule Curve and the USACE Rule Curve values for the same month. This figure shows the values for the month of march of Tayucai reservoir, Case 1, in which the values of the 2-Parameter Unrestricted A10 Rule Curve range between 352.3 and 380.0 masl, depending on the initial level in the reservoir; while the USACE Rule Curve is a constant at 370.4 masl. The range of levels between the two Rule Curves (shaded region) can be seen as a complementary solution space for a Rule Curve, which can be explored in diverse ways, either at random or methodically.

One form of exploration of the complementary solution space is by making a linear combination of level/storage from the two Rule Curves with the appropriate weights, plus the extreme values (maximum and minimum) of the range. Thus is it formed the Combined 2-Parameter-USACE Rule Curve, which can be tested by simulation to select the most appropriate curve for a reservoir and power-plant operating in a wide range of rated capacity.

Extreme Values

This combined extreme curve is represented by any of the following equations: Maximum Value: L _(C)=max(L _(2-P) ,L _(USACE))  a) Minimum Value: L _(C)=min(L _(2-P) ,L _(USACE))  b)

In which:

-   -   L_(C): Combined Rule Curve Level     -   L_(2-P): Level of the 2-Parameter Rule Curve     -   L_(USACE): Level of the USACE Rule Curve

FIG. 16 shows a graph of the combined maximum and minimum values for this example.

Linear Combinations

The linear combination of the level/storage is represented by any the following equations: Linear Combination on Level: L _(C) =w L _(2-P)+(1−w)L _(USACE)  a) Linear Combination on Storage: S _(C) =w S _(2-P)+(1−w)S _(USACE).  b)

The Combined Rule Curve level, L_(C), is then interpolated from the elevation-storage curve.

In which:

-   -   L_(C): Combined Rule Curve Level     -   L_(2-P): 2-Parameter Rule Curve Level     -   L_(USACE): USACE Rule Curve Level     -   S_(C): Combined Rule Curve Storage     -   S_(2-P): 2-Parameter Rule Curve Storage     -   S_(USACE): USACE Rule Curve Storage     -   w: weight assigned to the 2-parameter rule curve in the linear         combination, lies within the normalized range 0 to 1. The weight         can be the same throughout the year, or can vary with the         month-seasons

FIG. 17 shows several linear combinations on level and storage with weights 0.25, 0.50, and 0.75, as well as the original A10 and USACE Rule Curves. A linear combination on level is different than a linear combination on storage due to the nonlinearity of the level-area-storage curve of the reservoirs.

FIG. 18 shows the mean annual energy generated by the simulation of the combination of the extreme values (maximum and minimum) of the 2-Parameter Rule Curves (Unrestricted A 10 and Restricted A20) and the USACE Rule Curve. The maximum annual energy is generated by the use of the 2-Parameter Unrestricted A10 Rule Curve.

FIG. 19 shows the mean annual energy generated by the simulation of the liner combination on levels of the 2-Parameter Rule Curves (Unrestricted A10 and Restricted A20) and the USACE Rule Curve. The weights assigned to the levels of the 2-Parameter Rule Curves vary between 0.1 and 0.9, and are the same value throughout the year.

The simulated results show a smooth, non-linear variation between the energy generated by the 2-Parameter Rule Curves and the USACE Rule Curve as the weight assigned to the 2-Parameter Rule Curve varies from 0.90 down to 0.1. The maximum annual energy is generated by the use of the 2-Parameter Unrestricted A10 Rule Curve.

Linear Combination on Storage

FIG. 20 shows the mean annual energy generated by the simulation of the linear combination on storage of the 2-Parameter Rule Curves (Unrestricted A10 and Restricted A20) and the USACE Rule Curve. The weights assigned to the levels of the 2-Parameter Rule Curves vary between 0.1 and 0.9, and are the same value throughout the year.

The simulated results show a smooth, non-linear variation between the energy generated by the 2-Parameter Rule Curves and the USACE Rule Curve as the weight assigned to the 2-Parameter Rule Curve varies from 0.90 down to 0.1. The maximum annual energy is generated by the use of the 2-Parameter Unrestricted A10 Rule Curve.

In conclusion, the evaluation of the combination of the 2-Parameter and USACE Rule Curves indicates that the 2-Parameter Unrestricted A10 Rule Curve generates the maximum mean annual energy output.

Use of the 2-Parameter Rule Curve

The use of the 2-Parameter Rule Curve can increase total energy output from a single reservoir over the energy produced by using the USACE Rule Curve, while keeping the Firm Energy Yield of the reservoir, over a wide range of rated capacities in the power plant. The example is again Tayucai reservoir, in which the rated capacity of the power plant varies between 2,000 MW and 6,000 MW, with normal power pool level set at 380 masl and minimum power pool level set at 340 masl.

Characteristics of the Power Plant

FIG. 21 shows the general characteristics of the Power Plant, with rated capacities ranging from 2,000 MW up to 6,000 MW. The table contains the design point (head, flow, efficiency) and the range of net head (maximum and minimum). In this example the head varies with the rated capacity, and is calculated considering the level of the reservoir at the center of mass of the useful storage.

Energy Produced

USACE Rule Curve

FIG. 22 shows the energy output generated by Tayucai power plant using the USACE Rule Curve. It shows the continuous power, the firm energy yield, the secondary energy and average energy, as well as the plant capacity factor, ranging from 80% down to 25%. The maximum energy produced is around 15,240 GW-h/year with the installation of 2,480 MW in the power plant, at a plant capacity factor of 70%.

2-Parameter Unrestricted A10

FIG. 23 shows the energy output generated by Tayucai power plant using the 2-Parameter Unrestricted A10 Rule Curve, by itself or in combination with the USACE Rule Curve. It shows the continuous power, the firm energy yield, the secondary energy and average energy, as well as the plant capacity factor, ranging from 80% down to 30%. The maximum energy produced is around 17,460 GW-h/year with the installation of 4,600 MW in the power plant, at a Plant Capacity Factor of 43%.

The continuous power and firm energy generation does not change, however the average energy output increases by 2,240 GW-h/year which represents an increase of 15% over the maximum average energy and of 26% over the maximum secondary energy generated by the USACE Rule Curve alone.

2-Parameter Restricted A20

FIG. 24 shows the energy output generated by Tayucai power plant using the 2-Parameter Restricted A20 Rule Curve. It shows the continuous power, the firm energy yield, the secondary energy and average energy, as well as the plant capacity factor, ranging from 80% down to 30%. The maximum energy produced is around 16,980 GW-h/year with the installation of 4,400 MW in the power plant, at a Plant Capacity Factor of 44%.

The continuous power and firm energy generation does not change, however the average energy output increases by 1,740 GW-h/year which represents an increase of 11% over the maximum average energy and of 20% over the maximum secondary energy generated by the USACE Rule Curve alone.

Comparison

FIG. 25 shows the total energy (firm, average and secondary) generated by Tayucai Reservoir over the entire range of rated capacities (2,000 to 6,000 MW) analyzed. It shows the energy generated by the three rule curves: the USACE, and the envelope (maximum values) of the energy generated by the Unrestricted A10 and Restricted A20 in their diverse combinations (by themselves, as well as the maximum-minimum and the linear combinations with the USACE Rule Curve). The data is plotted in FIG. 26 (Mean Annual Energy) and FIG. 27 (Secondary Energy) over the entire range of rated capacity. Conclusion:

As can be seen, the Unrestricted A10 Rule Curve generates more energy from the reservoir for every rated capacity installed in the Tayucai Power Plant, from a low 2,000 MW (Plant Capacity 80%) to a high of 6,000 MW (plant capacity 30%).

The mean annual energy generated with the 2-Parameter A10 and A20 Rule Curves reaches a maximum value and remains fairly constant near the maximum value generated notwithstanding the increase in rated capacity of the power plant. 

The invention claimed is:
 1. A method for operating a hydropowered reservoir, the method providing an increase in average annual energy produced by the reservoir in comparison with energy generated by the reservoir when operated in accordance with a Rule Curve generated in accordance with a United States Army Corps of Engineers (“USACE”) operational strategy, without reducing annual firm energy yield for the reservoir, the method comprising: estimating continuous power and firm energy yield of the reservoir in accordance with the USACE operational strategy; determining a first Rule Curve for the reservoir in accordance with the USACE operational strategy; determining two variants of a composite monthly maximum parameter for generation of continuous power and firm energy at the reservoir; generating a second Rule Curve for the reservoir which considers the water level of the reservoir at the beginning of each month and one of the two variants of the composite monthly maximum parameter or a linear combination of the two variants of the composite monthly maximum parameter; and operating the reservoir in accordance with the second Rule Curve.
 2. The method of claim 1 wherein the composite monthly maximum parameter represents a monthly maximum effective potential energy withdrawn from, or stored into, the reservoir.
 3. The method of claim 1 wherein the step of generating a second Rule Curve takes into account one or more of the following for a given month: (a) expected rainfall into the reservoir; and (b) expected evaporation of water from the reservoir. 